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Abstract 



q ■ Application of Karhunen-Loeve decomposition (KLD, or singular value decomposi- 

tion) is presented for analysis of the spatio-temporal dynamics of wide-aperture vertical 
cavity surface emitting laser (VCSEL), considered as a thin-layer system. KLD tech- 
nique enables to extract a set of dominant components from complex dynamics of 
system under study and separate them from noise and inessential underlying dynam- 
ical behavior. Properties of KLD spectrum and structure of its main components are 
studied for different regimes of VCSEL. Along with the analysis of VCSEL, a brief 
survey of KLD method and its usage for theoretical and experimental description of 
nonlinear dynamical systems is presented. 

>- ■ Key words: thin films, spatio-temporal dynamics, VCSEL, singular value decomposition, 
^ ! Karhunen-Loeve decomposition 

O 

O ■ PACS numbers: 42.55.Px; 42.60.Jf 



UDC: 621.373.826.038+539.2 



O 

O . 

'55 ■ 1 Introduction 

H 

p ' Investigation of interaction of thin-film systems with laser radiation becomes quite topical 
during the last decade. This is mainly stimulated by advances of semiconductor technology, 



which enable to obtain multi-layer semiconductor structures with thickness close to or even 
less than wavelength of visible radiation. In particular, wide-aperture vertical cavity surface 
& ', emitting lasers (VCSELs) are wide-spread optical systems of communication and information 
processing. Complex spatio-temporal regimes of VSCEL operation, in principle, opens new 
possibilities for information processing (for example, "chaotic" encoding for secure applica- 
tions). On the other hand, it is far from full understanding, what mechanisms form spatial 
and temporal structure of radiation in VCSELs. 
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In the paper, we investigate the spatiotemporal dynamics of VCSEL using the system of 
differential equations describing the dynamics of broad area VCSEL [1]. The equations are 
derived using spin flip model [2,3] and take into account polarization of light and complex 
cavity of VCSEL, including Bragg reflectors. The propagation of radiation through the 
VCSEL cavity is calculated in an approximation of thin nonlinear layer, which allows to 
dismiss the diffraction of light in the active medium. Calculations of the current intensity 
J near threshold and for some values far from threshold are performed for the case of zero 
phase anisotropy, low amplitude anisotropy and homogeneous spatial profiles of refraction 
index and injection current. 

Obtained sets of spatiotemporal data are analyzed using Karhunen-Loeve decomposition 
(KLD, also known by several other names [4]). Such technique is introduced in the beginning 
of 20th century for description of random functions. Then this method have found numerous 
applications in such areas as pattern recognition, turbulence, meteorology, coherence theory 
etc. Optimal properties of KLD enables to extract only a few main components from the 
whole complex dynamics of system under study. 

The rest of article is organized as follows: in the next section we outline the method 
of or Karhunen-Loeve decomposition and its main characteristics. In the third section we 
provide the mathematical model of VCSEL is presented. In the following section we analyze 
of complex dynamics of VCSEL by KLD method. In particular, change of the decomposition 
spectrum for the values of the current density J near lasing threshold and for some values 
above threshold is discussed. The last section contains a conclusion and outline of future 
tasks. 



2 Karhunen-Loeve decomposition 

In its simplest form, suitable for our purposes, the KLD method is formulated as follows: 
given some (in general, complex) field of two variables u(r,t) with r = (x,y), one tries to 
find its decomposition onto purely temporal and spatial modes: 

oo 

u(r, t) = ^ A, flj(t) <f>i(r). (1) 
i=i 

with two orthonormality conditions 

/ dv <Pi(r) 0*(r) = I dt Oi{t) a*(t) = % (2) 

Js Jt 

where T is a time interval and S is an area on which we want to analyze field. This field is 
considered as known, e.g. from experiment, or from some kind of model - either analytical 
or numerical. 

The physical sense of representation is extraction of spatial distributions which os- 
cillate in time as a whole. The values |Aj| 2 give the part of 'energy' carried by i-ih mode 
in average. Physical origin of u(r,t) is of little importance — it could be electrical field or 
intensity of some kind of radiation, velocity profile of flow or even simply set of pictures u(r) 
numbered by second variable t [5,6]. In our case, u(r) is the slowly varied complex envelope 
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of the optical field. One of the main parameters of presentation {IJ is the number of terms 
N e , such that their sum contains not less that some prescribed part of the whole energy of 
u(r,t) 

N e oo 

£|A*| 2 >(l-e)X>| 2 - (3) 

i=l i=l 

The other components with i > N e bears noise or unimportant dynamics and so could be 
excluded from consideration (on given interval of time and spatial domain). 

As it could be shown, the decomposition functions could be found from two eigenproblems 
for integral equations 

|A,| 2 ^(r)= / drtfaWKsM (4) 
Js 

and 

\k\ 2 a l (t) = J dt 1 ai(t')K t (t,t') (5) 
where kernels are correlation function 

#s(r,r') = f dtu(r,t)u(r',t) (6) 
Jt 

and some kind of temporal correlation function, averaged over space 

K s (r,r') = ^dtu(r,t)u(r',t). (7) 

Equation (J3J) with kernel (J7J) is usually referred to as "method of snapshots" of "method of 
strobes" [4]. 

On the other hand, decomposition of type ((H), (j2J corresponds to singular- value decom- 
position [7] of "matrix" u(r,t). As far as experimental (or numerical calculation) data is 
always is a kind of matrix of numbers. Appropriate discrete decomposition may be effectively 
calculated using standard svd routine available in umber of mathematical packages. 

It should be noted, that whole decomposition could be found from just one eigenproblem, 
while the dual basis, temporal or spatial, is found from projection 

\ai(t) = dr 0*(r) w(r,t), 
Js 

Xiipi(r)= / dt a*(t)u(r,t). 
Jt 

Use of eigenfunctions related to the investigated field cause main positive sides of 
Karhunen-Loeve expansion: among all sums of type (0) with finite number of terms, the 
representation in terms of eigenfunctions ensures minimal least square error and the max- 
imal capture of "energy." The KLD method proves its power on number of problems in 
different areas of physics and other sciences. On the other hand, its main drawback again 
related with use of eigenfunctions: analytical solution of equations P|. is known only for 
very few special cases, the numerical solution often require too much resources and is inca- 
pable to provide all the information about system dynamics, especially near critical points. 
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To this point, it is important to look for methods of analysis, which provide information 
about spectrum of eigenvalues |A| 2 without calculation of decomposition itself [8]. 

However, calculation of Karhunen-Loeve decomposition provides valuable information 
about the details of complex process. Solution of eigenvalue problem is usually much more 
easy task than to study other, "standard" parameters of chaotic systems, such as Lyapunov 
exponents or fractal dimensions. In most cases singular-value analysis enables to select just 
the few most important components from whole spatio-temporal dynamics and to study 
their behaviour. 
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3 Short description of the VCSEL model 

Response of the active medium to the radiation in VCSEL is in a semiclassical approximation 
is described by the spin-flip model [1-3], taking into account vector character of field: 

is] P - SzAE, 
/ 3h 

- -4 (E*P - EP*) , (8) 

2h 

L (e*P' - EP'*) . 

where P is the polarization of the two level centers, P' is the vector with components 
(P y , —P x )] N is the total population difference between the conduction and the valence bands 
and Nq is its transparency value; n is the difference between the population differences for 
the two allowed transitions between magnetic sublevels, 5 = u g — oj c is the detuning between 
the bandgap frequency u) g and the cavity resonance uj c \ ^ s is the decay rate between the 
magnetic sublevels; T\ and T 2 are the relaxation times for the total population difference 
and the polarization correspondingly; \d\ is the absolute value of the dipole momentum of 
the transition (we suppose it is the same for both transitions); J is the pump parameter, 
and 

A= (N-N in \ 

\ -in N — N J ' {J> 

We will use equations (JHj) with adiabatically eliminated polarization. The procedure of 
such adiabatic elimination is described in details in [1] and allows to take into account the 
asymmetry of the gain line using so called linewidth enhancement factor a. In addition, this 
procedure allows to avoid the short-frequency instabilities intrinsic to the straightforward 
adiabatic elimination procedure for the spatially extended lasers. In this approximation the 
polarization of the active medium is defined as follows: 



o ^(i-a)D£E, (10) 

here D = N — N , the operator £ ~ £(k±) = 1/[1 + T|(<5 — Q(k±)) 2 }, where Q(k±) is a cavity 
resonance frequency for the tilted wave with definite k±, and the tilde means an equivalence 
in the sense of transverse Fourier transform. 
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The propagation of radiation through the VCSEL cavity is calculated in an approximation 
of thin film active medium, which allows to neglect the diffraction in an active layer [1]. It 
gives us the following relation: 



where Ei is the field incident into the active medium, E t is the field outgoing from the active 
medium. Operator F = p exp(2ikL + i(A±/k)L) T R, where p describes absorption in the 
linear medium between the active layer and reflector, L is the thickness of the spacer layer, 
is the Laplasian in the transverse plane, 



is an polarization anisotropy matrix, with 7 a , 7 P being the amplitude and phase anisotropy 
parameters correspondingly. Operator R describes the reflection of plane waves from the 
Brag reflectors [1]. 

For numerical simulations we used the following parameters: 7 S = 100, 7 a = 0.1, 7 P = 0.0, 
5 = 0.006, a = 3.0. For chosen parameters the lasing threshold is Jo — 0.730 which is needed 
for the following discussion (see [1] for the detailed description of threshold conditions). 

4 Analysis of VCSEL dynamics 

We consider first the simplest homogeneous case with periodic boundary conditions and with 
small amplitude anisotropy. The resulting spatiotemporal regime near threshold (J = 0.740) 
is regular and consists of stripes that weakly oscillate near certain equilibrium state (the 
directions of oscillations is shown in Fig^J) by arrows). The contrast of stripes also changes 
during the evolution. The KL-spectra of eigenvalues Xj of the spatiotemporal regime for this 
case is presented in Fig^i. It is evident, that only two modes are the most significant, while 
all the others could be safely treated as zeros. 

The spatial KL- modes appear to be stripes (FigEJ). Their maxima don't coincide and 
thus a spatial phase shift leads to orthogonality and to spatiotemporal dynamics. Hence, 
the averaging over a large period of time leads to stripes with smaller contrast. Time depen- 
dencies of the KL-modes (see Fig. EJ) also characterized by regular (oscillating) dynamics, 
except only some transient stage. 

However, a slight enlarge of the injection current (up to J = 0.800) leads to chaotic time 
dependency the field. Nevertheless, in this case the KL-spectra has the same form as one 
near threshold Fig. and, moreover, the spatial KL-modes are the same as in Fig. El too. 
However, their temporal behavior is sufficiently different (see FigEJ). 

Therefore, the chaotization of the regular dynamics appears due to mechanism, which 
is not connected with excitation of the long-wavelength inhomogenities, either Eckhaus or 
zig-zag type (i.e. in x or y directions). 

Further increase of the injection current (J = 1.000) leads to more disordered dynamics. 
More than two active modes are present in in KL-spectrum (Fig. EJ). The spatial structure is 
also changing in this case, and modulation in x direction now is accompanied by modulation 
in y direction (see Fig. EJ). 



Ej — F E t . 



(11) 




(12) 



5 



5 Conclusion 



In summary, our calculations have shown, that during transition of VCSEL's from the regular 
behavior to spatio-temporal chaos, it still can be described by superposition of just a few 
modes with relatively simple structure (for moderate values of the injection current). Increase 
of order parameter (injection current J) leads to activation of some new modes, with new 
features of transversal and temporal dependence. 

Hence, the observed chaos is not truly "spatio-temporal". Complex dynamics in time 
domain (both for the electomagnetic field itself and for KL-modes) is accompanied by just 
very simple spatial structure of modes. Moreover, the whole dynamics is described by only a 
few components. This fact, together with importance of VCSELs in modern optical commu- 
nication, enables to suppose, that chaotic regimes could be effectively controlled by adjusting 
parameters of a system. 
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Figure 1: (a) - The KL spectra for J=0.740. (b) - The zoomed snapshoot of the dynamics. 
The arrows shows the direction of oscillations(the amplitude of oscillations is smaller then 
the length of arrows). 
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Figure 2: The first 2 spatial KL-modes for J = 0.740. 
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Figure 3: The temporal evolution of KL-modes shown in Fig. |21 
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Figure 4: The first 2 temporal KL-modes for J=0.800. 
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Figure 5: The KL spectra for J=1.000. 
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Figure 6: The first 6 spatial KL-modes for J=1.000. 



